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In this Letter, we present a theoretical analysis of the acoustic transmission through two-dimensional 
arrays of straight rigid cylinders placed parallelly in the air. Both periodic and completely random 
arrangements of the cylinders are considered. The results for the sound attenuation through the 
periodic arrays are shown to be in a remarkable agreement with the reported experimental data. 
As the arrangement of the cylinders is randomized, the transmission is significantly reduced for a 
wider range of frequencies. For the periodic arrays, the acoustic band structures are computed by 
the plane-wave expansion method and are also shown to agree with previous results. 
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When propagating through media containing many 
scatterers, waves will be scattered by each scatterer. The 
scattered waves will be scattered again by other scatter- 
ers. This process is repeated to establish an infinite recur- 
sive pattern of rescattering between scatterers, forming 
a multiple scattering process 0,^]. Multiple scattering 
of waves is responsible for a wide range of fascinating 
phenomena, including twinkling light in the evening sky, 
modulation of ocean ambient sound || , acoustic scintilla- 
tion from fish schools [0]. On smaller scales, phenomena 
such as white paint, random laser [j5j, electron transport 
in impured solids |6| are also results of multiple scatter- 
ing. When waves propagate through media with periodic 
structures, the multiple scattering leads to the ubiquitous 
phenomenon of band structures. That is, waves can prop- 
agate in certain frequency ranges and follow a dispersion 
relation, while within other frequency regimes wave prop- 
agation is stopped. The former ranges are called allowed 
bands and the latter the forbidden bands. In certain sit- 
uations, the inhibition of wave propagation occurs for all 
directions, leading to the phenomenon of complete band 
gaps. 

The wave dispersion bands are first studied for elec- 
tronic waves in solids, providing the basis for understand- 
ing the properties of conductors, semi-conductors, and 
insulators uj. In late 1980s, it became known that such 
a wave band phenomenon is also possible for electro- 
magnetic waves in media with periodically modulated 
refractive- indices jsj . Since then, optical wave bands have 
been extensively studied, yielding a rich body of litera- 
ture. The theoretical calculations have proven to match 
well with the experimental observations . 

In contrast, research on acoustic wave band struc- 
tures has just started. Although theoretical computa- 
tions of band structures have been documented for peri- 
odic acoustic structures ]To| , the experimental work was 
only recent, and to date only a limited number of mea- 
surements has been reported. One of the first observa- 
tions was made on acoustic attenuation by a sculpture 
pT|. The authors obtained a sound attenuation spec- 



trum, which was later verified by the band structure com- 
putation |lj]. Recently, acoustic band structures have 
been further measured for acoustic transmission through 
two-dimensional (2D) periodic arrays of rigid cylinders 
placed in the air The authors demonstrated the 

properties of sound attenuation along two high-symmetry 
directions of the Brillouin zone of the arrays. They also 
observed a peculiar effect of deaf bands; within the bands, 
in spite of non-zero band states, wave propagation is pro- 
hibited due to particular symmetry of the states |l3| . 

The main purpose of this Letter is to provide a theo- 
retical investigation of sound transmission by 2D arrays 
of rigid cylinders in air in line with the experiment of 
p3[ , thus providing a direct comparison of attenuation 
between theory and experiment. Such a direct compar- 
ison is relatively scarce in the literature. We note that 
the comparison between the attenuation spectrum and 
the dispersion bands is of indirect nature, as the two are 
not necessarily one to one correspondent, as seen, for in- 
stance, when some seemingly allowed bands are actually 



deaf to wave transmission |13|. This will be further clari- 
fied in the later results. Another goal is to study how the 
randomness affects the acoustic transmission, so to gain 
some insight into the connection between the forbidden 
bands and wave localization in 2D fl4j| . For the purposes, 
we adopt a self-consistent multiple scattering theory [15l . 

Consider N straight cylinders located at with i = 
1, 2, • • • , N to form either a regular lattice or a random 
array perpendicular to the x — y plane; the regular ar- 
rangement can be adjusted to comply with the experi- 
ment [p"3f . The cylinders are along the z-axis. An acous- 
tic source transmitting monochromatic waves is placed 
at f s , some distance from the array. The scattered wave 
from each cylinder is a response to the total incident 
wave composed of the direct wave from the source and 
the multiply scattered waves from other cylinders. The 
final wave reaches a receiver located at f r is the sum of 
the direct wave from the source and the scattered waves 
from all the cylinders. Such a scattering problem can 
be formulated exactly in the cylindrical coordinates, fol- 
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lowing Twersky While the details are in |T(|, the 

essential procedure is presented below. 

The scattered wave from the j-th cylinder (j = 
1, 2, • • • , N) can be written as 

oo 

p s {r,f j ) = ™A^HV>{k\? (1) 

n— — oo 

where i = a/ - 1) Hn is the n-th order Hankel function of 
the first kind, (pp-p, is the azimuthal angle of the vector 
f — fi relative to the positive x-axis. The total incident 
wave around the i-th cylinder (i = 1, 2, • • • , N; i ^ j) is 

N 

p\nif) = Poif) + E P'(f^i)> ( 2 ) 
i=i ,3& 

which can be expressed again in terms of a modal series 

oo 

Pln(f)= E BiJ n (k\f-fl\)e in ^-n. (3) 

n— — oo 

The expansion is in terms of Bessel functions of the first 
kind J n to ensure that p\ n {f) does not blow up as r — > i\. 

To solve for A % n and B l n , we express the scattered wave 
Ps{t, fj), for each j ^ i, in terms of the modes with 
respect to the i-th scatterer by the addition theorem for 
the Bessel functions |l7]]. The resulting formula for the 
scattered wave p s (r , fj ) is 

oo 

Ps(r,fj)= E C^J n (k\r-rl\)^-^, (4) 

71— — OO 

with 

oo 

C#*= E mAjH&m-r- , j \V < - l ~ n) *' i -*>- (5) 

l — — oo 

The direct incident wave around the location of the i-th 
cylinder can be expressed in the Bessel function expan- 
sion with respect to the coordinates centered at fi 

oo 

Po(f)= E S\mf-fl\)e u ^\ (6) 

I — — OO 

with the known coefficients 

Si = mH^(k\fl-f a \)e- il M. 

Matching the coefficients in equation (||), using equa- 
tions (||), (Q) and (||), we have 

N 

B i n =si+ e c n- ( 7 ) 

At this stage, both the 5^ are known, but both B l n and 
are unknown. Boundary conditions will give another 



equation relating them. The boundary conditions state 
that the pressure and the normal velocity be continuous 
across the interface between a scatterer and the surround- 
ing medium. After a deduction, we obtain 

4 = mTiAi, (8) 

where 

g i h i J^(ka i )J n (ka i /h i ) - J n (fca') J^kat/h*) 

(9) 

Here the primes refer to taking derivative, a 1 is the radius 
of the i-th cylinder, g l = p\/p is the density ratio, and 
h L = k/k\ = c\/c is the sound speed ratio for the i-th 
cylinder. 

The unknown coefficients A % n and B° n can be inverted 
from Eqs. (§), (0), and §). Once A l n are determined, the 
transmitted wave at any spatial point is given by 

JV oo 

p(r)=Po(r)+E E ^KH^(k\r-n\)e in ^. 

i — 1 n— — oo 

(10) 

The acoustic intensity is represented by the squared mod- 
ule of the transmitted wave. When the cylinders are 
placed regularly, we can also obtain the band structures 
by the plane wave method well documented in JhJ . The 
programs used are identical to that used for computing 
the acoustic bands in the regular arrays of air-cylinders 
in water Jl8| 

Numerical computation has been carried out for the 
experimental situations Jl3[ and also for an random ar- 
ray of the cylinders. In the simulation, all the cylinders 
are assumed to be the same, in accordance with the ex- 
periment. Moreover, the radii of the cylinders and the 
lattice constants are also taken from the experiment. Sev- 
eral values for the acoustic contrasts between the cylinder 
and the air were used in the initial stage of computation. 
We found that the results are in fact insensitive to this 
factor as long as the contrasts exceed a certain value. 
This agrees with the experimental observation. In the 
results shown below, we use g = h = 20 as the values for 
the acoustic contrasts. In the computation, we allow the 
number of the total cylinders to vary from 100 to 500, in 
line with the experiment. In the particular results shown 
later, we assume that the cylinders arc placed within a 
rectangular area of 8 x 40 of lattice domain. The source 
and receiver are placed about one lattice constant away 
from the long side of the array so to minimize the effect 
due to the finite sample size. As we do not know the spec- 
ifications for the transmitter and the receiver, we assume 
an omni-directional transmitter as the acoustic source lo- 
cated on one side of the array of the cylinders, whereas 
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an omni-directional receiver is placed on the other side 
to receive the propagated waves. 

In Fig. 1, we show the relative attenuation (~ — In \p\ 2 ) 
spectra for various square lattices for acoustic trans- 
mission along the [100] direction. The parameters for 
the four cases considered are adopted from the experi- 
ment 13 1 . We observe a robust attenuation peak located 
around 1.5 kHz for all the situation. By eye inspection 
of Fig. 1 and Fig 1 in the agreement between the 
theoretical and the measured results are good, particu- 
larly in view of the finite dimension of the arrays and 
no adjustable parameters. The height of the attenuation 
peaks depend on the locations of either the receiver and 
the transmitter, on the outer boundary of the cylinder ar- 
rays, and on the number of cylinders. Nevertheless, the 
theoretical results describe quantitatively well the obser- 
vation. A slight difference appears, however, for the case 
with cylinders of diameter 1 cm at the filling of 0.006. In 
our results, a small attenuation peak occurs at about 1.5 
kHz, but is absent from the experiment. This discrepancy 
may be attributed to a couple of reasons: the theoretical 
setting does not match exactly that in the experiment 
and perhaps the attenuation peak is too small to be ob- 
servable. In Fig. 1 we also observe some peaks located at 
higher frequencies as observed in the experiment. These 
peaks are sensitive to the arrangement of the transmitter 
and receiver, and the number of the cylinders. They are 
not because the frequencies are within a stop band. 
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1. Theoretical results for the relative acoustic atten- 



uation along [100] as a function of frequency for various ar- 
rangements of cylinders. The parameters for the cylinder ar- 
rays are taken from the experiment. 



Fig. 2 shows the attenuation (relative) spectra for the 
case of cylinders of diameter 3 cm in a square lattice with 
lattice constant 11 cm. The corresponding band struc- 
ture is also depicted. The attenuation along the [100] and 
the [110] directions are represented in solid and dotted 
lines on the right panel respectively. Comparing Fig. 2 



with the experimental Fig. 2 in [[L3|, we see that the at- 
tenuation peak along [100] coincides almost exactly with 
the experimental data in the frequency range between 
1.38 and 1.70 kHz. In this particular case, the attenua- 
tion peaks are also roughly in the same order of magni- 
tude as the observation. Along the [110] direction, the 
attenuation due to the deaf bands are also nicely recov- 
ered by the theory. That is, the two bands predicted 
by the band structure computation (the second and the 
third bands in Fig. 2 for the range from M to L) are 
actually deaf and wave propagation is prohibited within 
these two bands. Similar results have also been reported 
for 2D photonic band gap materials [p|Jl9|]. Further sim- 
ulation shows that though the height of the attenuation 
peaks may vary as the outer boundary of the cylinder ar- 
rays changes, the overall shape of the attenuation spectra 
remain unchanged. 




Band Structure 



Sound attenuation (dB) 



FIG. 2. Right panel: Relative acoustic attenuation vs. 
frequency for a square array of cylinders with periodicity 11 
cm. The radius of the cylinders is 3 cm. Left panel: The band 
structures computed by the plane wave expansion method. 



We also notice that the width of the attenuation peak 
along [110] is a little wider than the observation. In ad- 
dition to the aforementioned reasons, this discrepancy 
could be due to the fact that in the present simulation, 
the cylinders are assumed to be in the open air, while the 
experiment was performed in a chamber which may some- 
what still reflect sound. Furthermore, the exact number 
and the setting of the cylinders in the experiment are 
also not known from the literature. Other contributions 
to the discrepancy may result from the different acoustic 
source and receiver used in the theory and experiment. In 
spite of these limitations, the match between the theory 
and experiment is quite encouraging. 

The band structure shown in Fig. 2 is obtained by us- 
ing the usual plane- wave method pf]| . Here it is shown 
to agree nicely with the band structure obtained by the 
variation method calculation. With the parameters in 
Fig. 2, wave propagation in different directions is in- 
hibited within different frequency regimes. No complete 
band gap is observed. The overlap of attenuation peaks 
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along different directions, an indication of the complete 
band gap, can observed when cylinder filling factor ex- 
ceeds certain values [|f3| ; in this case, the matrix inversion 
in the multiple scattering computation becomes costly. 




Band Structure Sound attenuation (dB) 

FIG. 3. The right panel shows the relative acoustic atten- 
uation vs. frequency. Here the comparison is made between 
the results from a square array with periodicity 11 cm (solid 
line) and from a complete random array of cylinders (dot- 
ted line). The radius of the cylinders is 3 cm. Left panel: 
The band structures computed by the plane wave expansion 
method. 



Now consider the effect of randomness on the acoustic 
transmission. Here we take the case described by Fig. 2 
as the example. While keeping the same cylinder fill- 
ing and the outer boundary of the array, we allow the 
cylinders to be distributed completely randomly within 
the area. The relative attenuation is computed again. 
The results are shown in Fig. 3, with comparison to the 
results for wave propagation along [100] for the corre- 
sponding square latter array. The following features are 
evident. At low frequencies, the introduced disorder does 
not affect the transmission for the given size of the sam- 
ple. The randomness reduces the transmission basically 
for all frequencies above a critical frequency of 700 Hz, 
except for the frequencies within the lowest stop band 
along [100] within which though the transmission is still 
inhibited, the disorder actually reduces the attenuation 
purely due to the stop-band effect; such a feature is also 
observed in other acoustic systems po| . As we increase 
the sample size, the critical frequency tends to become 
smaller, implying larger ranges of inhibition. Note that 
in order to compute the transmission accurately at lower 
frequencies, larger sample sizes are required and the com- 
putation would become costly. The result of the severe 
reduction in transmission for a wide range of frequencies 
is remarkable and has significant relevance to the fun- 
damental problem of Anderson localization, a concept 
originally introduced to explain the conductor-insulator 
transition induced by disorders in electronic systems , 
because such a reduction is a precursor to the localiza- 
tion phenomenon. The results also imply that random 



arrays of rigid cylinders are good candidates in filtering 
audible noise. 

In summary, here we have presented a theory for acous- 
tic transmission through arrays of rigid cylinders in the 
air. The theory are applied to the experimental situ- 
ations of regular arrays, yielding favorable agreements. 
The theoretical results verify the existence of the deaf 
bands. The results are subsequently extended to the 
case of random cylinder arrays. We found that wave 
propagation is significantly reduced by randomness for a 
wide range of frequencies. This feature makes the rigid 
cylinders in the air an ideal system for theoretical and 
experimental studies of wave localization. 
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